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Abstract 

Set differential equations are usually formulated in terms of the 
Hukuhara differential, which implies heavy restrictions for the nature 
of a solution. We propose to reformulate set differential equations as 
ordinary differential equations in a Banach space by identifying the 
convex and compact subsets of 1R4 with their support functions. Us¬ 
ing this representation, we demonstrate how existence and uniqueness 
results can be applied to set differential equations. We provide a sim¬ 
ple example, which can be treated in support function representation, 
but not in the Hukuhara setting. 
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1 Introduction 

A set differential equation is an equation of the form 

D H A(t) = f(t,A(t)), A(0) = A Ol (1) 
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where t i->- A{t) is a curve in the space /C c (R d ) of nonempty convex and 
compact subsets of R cZ , the right-hand side is a mapping 

/ : [0,T] x /C c (R d ) —> /C c (R d ), 

and D H A(t) is the so-called Hukuhara differential of the curve at t e (0,T). 
Set differential equations have been investigated in a considerable number 
of papers. For an overview of the literature we refer to m- The usage of 
the Hukuhara differential in (JTJ) implies heavy restrictions for the nature of 
a solution, which can, e.g., only grow in diameter, but not shrink, see m 
Proposition 1.6.1]. 

Recently, there have been attempts to modify the underlying Hukuhara 
difference with the aim to allow for a more flexible behavior of solution curves, 
see mm and the references therein. The resulting differential is called the 
second type Hukuhara differential. In this setting, solution curves of ([IJ) can 
shrink, but not grow. There exist, however, curves in /C c (R d ) with d > 2, 
which expand in some space directions and contract in others simultaneously. 
Both Hukuhara-based approaches fail to capture this behavior. 

The Hukuhara differential is not the only approach to handle set evolu¬ 
tions. In particular, we would like to mention an abstract framework named 
Mutational Analysis , which has been presented in [T] and further developed in 
[12] • It generalizes evolution equations from vector spaces to metric spaces 
and can not only handle evolutions in /C c (R d ), but also in spaces of more 
general sets such as the compact subsets of R d . 

The aim of the present paper is to show that a large family of evolutions 
in /C c (R d ), containing the problems investigated in [IT] and [121 E], can be 
written and treated as ordinary differential equations in a Banach space with 
the usual Frechet derivative in time. We do not apply the apparatus from 
P and 112], but obtain very satisfactory results by exploiting the intrinsic 
features of of /C c ( R d ). 

Identifying convex sets with their support functions yields an embedding 
of the space /C c (R d ) into the Banach space of continuous real-valued 

functions on the sphere, see [9]. As it is well-known that any Hukuhara dif¬ 
ferentiable curve is Frechet differentiable in support function representation, 
see [31 Lemma 4.1], it seems natural to consider set differential equations in 
support function representation 

f t °A{t) = ch4(o)=<A4 0 , (2) 
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where t H > A(t) is a curve in /C c (R d ), t H y crA(t) is a curve in C(§ d_1 ), and 
is the Frechet differential of the curve at £ € (0, T). 

There are some technical difficulties when standard results on ordinary 
differential equations are applied to equations of type (j2]). As we have to 
guarantee that solutions stay in the manifold £ C (7(§ d_1 ) of all support 
functions associated with sets from /C c (R d ), we have to understand the struc¬ 
ture of the tangent cone T s ((t) to £ at any o G £. To transfer existence and 
uniqueness theorems for ordinary differential equations in Banach spaces with 
non-Lipschitz right-hand side to 02]) , we need compactness properties of £ 
and a characterization of the semi-inner product on (C(S d_1 ), || • ||oo). Some 
of these preliminary results can be taken from the literature, others are devel¬ 
oped in the present paper. In particular, we give a geometric interpretation 
of the one-sided Lipschitz condition in JC c (JR d ), which is a surprisingly mild 
condition on the behavior of /. 

The organization of the paper is as follows. In Section [2] we collect ba¬ 
sic definitions and the preliminary results mentioned above, which we use in 
Section E] to transfer standard existence and uniqueness results to (J2|). In Sec¬ 
tion [4] we briefly show that second-type Hukuhara differentiable curves are a 
special case of ([2]) . The example discussed in Section [5] illustrates that both 
Hukuhara approaches fail to capture very simple dynamics in /C c (R 2 ), while 
the support function calculus is applicable and yields reasonable solutions. 


2 Preliminaries 

After introducing basic notation in Section 12.11 we will collect some known 
results about support functions and tangent cones in Sections 12.21 and 12.31 
Section 12.41 investigates duality concepts, which are ingredients for standard 
results on ordinary differential equations in Banach spaces, in the particular 
case of set differential equations. 

2.1 Basic definitions 

Let Rq" be the set of all nonnegative real numbers. Throughout this paper, 
S^ 1 C R cZ will denote the sphere w.r.t. the Euclidean norm || • || : R d —> Rq", 
and the modulus will be denoted | • | : R —> Rg . Let C(S d_1 ) be the space of 
continuous real-valued functions on equipped with the maximum norm 
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II • lloo : C^" 1 ) -> Rj. If (X, || • || x) is a normed space, x E X and r > 0, 
then 

® r (x) := {x' E X : ||x' — x|| < r} 

is the closed ball of radius r centered at x. 

The nonempty compact subsets of R d will be denoted /C(R d ), and /C c (R d ) 
will stand for the nonempty convex and compact subsets of R d . For any 
A E R and A, B E /C(R d ), let 

A T B — {o T b \ cl E A, b E R} and XA ‘ — {An : <x E A } 
denote Minkowski addition and multiplication. For any A,B E /C c (R d ), let 
dist(7l, B ) := sup inf ||a — b||, 

aeA b ^ B 

dist #(A B ) := max{dist(A B ), dist(R, A)} 

denote the one-sided and the symmetric Hausdorff distance. For a,b E R d , 
we write dist(a,R) and dist(Al,6) instead of dist({a},R) and dist(7l, {6}). 
The projection of a point a E R d to a set B E /C(R d ) is the nonempty set 

proj B (a) \= {b E B : ||a — 6|| = dist(a, B)}. 

When B E /C c (R d ), then a i— y proj B (a) is a single-valued mapping, see [5] 
Lemma 7.3], and it follows from [5] Proposition 7.4] that this mapping is 
1-Lipschitz. 

We associate convex and compact subsets A E /C c (R d ) with their support 
functions 

a a ■ § d-1 —> R, cta(p) := sup(p, a). 

a£A 

Sometimes, it is useful to consider their positive homogeneous extensions 
a a : R cZ —> R, &a(p) := sup(p, a), 

aeA 

which obviously coincide with cta(-) on § d_1 . We dehne 

£(R d ) := : 71 G /C c (R d )} 

to be the set of all support functions of convex and compact subsets of R d , 
and we set 

S(R d ) := £(R d ) - £(R d ) = {a A ~ cr B : A, B E /C c (R d )}. 
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2.2 Elementary facts about support functions 

The following proposition is Corollary 13.2.2 from |l~5j . 

Proposition 1. A bounded function a : S^” 1 —>■ R is a support function of 
some A E /C c (R d ) if and only if its positive homogeneous extension d : R d —* 
R is convex. 

Recall that every convex function d : R d —> R is continuous (see [ 151 
Theorem 10.1]). We may therefore interpret the set E(R d ) of all support 
functions as a subset of C(§ d_1 ). 

The following facts are well-known (see IM). 

Proposition 2. If A, B E /C c (R d ) and A > 0, then 

a ) &a+b — & a + <Jb an d ct\a = Aa a, 

b) dist (A,B) = max peBl ( 0 ) (d A (p) ~ d B (p)), 

c) dist h (A,B) = maXpggd-i | a A (p) - a B (p) I- 

In particular, S(R d ) is a convex subcone of C^ -1 ). 

The cone S(R d ) is locally compact. 

Proposition 3. The cone £(R d ) is closed as a subset of C{and for 
any a E £(R d ) and r > 0, the intersection £(R d ) fl B r (cr) C is 

compact w.r.t. the maximum norm. 

Proof. Let a E C^S^ 1 ), and let (a n ) nG M C £(R d ) be a sequence of support 
functions with \\cr n — cr||oo —> 0 as n —* oo. By Proposition Q] the extensions 
d n are convex. Hence, we have for any A E [0,1] and x, y E R d that 

d(\x + (1 - A )y) ^d n (\x + (1 - A )y) 

<A d n (x) + (1 - X)d(y) -d A d(x) + (1 - X)d(y) 

as n —y oo, so that d is convex. Therefore, Prosition [T] implies that a E 
E(R d ). 

By Blaschke’s selection theorem, see m Chapter 4], the set S(R a! ) fl 
®IMIcx>+r-(0) is compact. As £(R d ) fl B r (cr) is the intersection of two closed 
sets, it is a closed subset of the compact set E(R d ) fl ®|| o -|| oo +r(0), and hence 
compact. □ 
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2.3 Tangent cones 

We are interested in C'(S d_1 )-valued solutions of differential equations that do 
not leave £(R d ). The concept of tangency is central for existence theorems 
under state constraints. 

Definition 4. Let X be a normed space, K C X a set and x G K . Then the 
tangent cone to K at x is given by 

T k (x) := {« 6 I : liminf h _1 dist(x + hv, K) = 0}. 

h \0 

The following proposition is Lemma 4.2.5 in [2]. ft will later be used to 
characterize tangency to the convex cone £(R d ). 

Proposition 5. If X is a normed space and K C X is a convex cone, then 
T k (x) = K + Ra; for all x G K. 

2.4 The semi-inner product for support functions 

In Section [21 we will apply a uniqueness theorem for ordinary differential 
equations in Banach spaces to set differential equations in support function 
representation. Its main ingredient is a one-sided Lipschitz condition, which 
is given in terms of a so-called semi-inner product. Therefore, we investigate 
in the present paragraph how this product acts on S(R d ) C C^” 1 ) and 
what this action means for the corresponding elements of /C c (R cZ ). 

Definition 6. For any Banach space X with dual space X*, the duality map 
J : X =4 X* is given by 

J(x) = {x* G X* : x*(x) = ||a;|LY = \\x*\\ 2 x *}. 

The mapping (•, •)_ :IxIgK defined by 

(x,y)~ = inf {y*(x) : y* G J{y)} 

is called a semi-inner product. 

Consider the Banach space X = C(M), where M is a compact metric 
space and C(M ) denotes the space of all continuous real-valued functions 
on M equipped with the maximum norm. Let B(M) denote space of all 
signed Borel measures on M, and let B(M) + denote space of all positive 
Borel measures on M. 
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Proposition 7 (Jordan decomposition). For any p G B(M), there exists a 
unique pair (pp,Pn) £ B + (M) x B + (M ) supported on Borel sets P, N C M 
such that p = pp — p^ and M is the disjoint union of P and N. 

For a proof, see Theorem 10 and Corollary 11 in Section III.4 of [7J. 

As a consequence, the total variation of a signed Borel measure is well- 
defined. 

Definition 8. The total variation of a Borel measure p G B(lR. d j with Jordan 
decomposition pp + pn = p with associated Borel sets P U N = M is defined 
by 

Ml p(P) + Hn(N). 

It is well-known that the dual space of (C(M), || • ||oo) is ( B{M ), |||-|||), 
which follows from the Riesz representation theorem, see Theorem IV. 6.3 in 

m 

We will now characterize the duality map on C(M). For a given function 
/ G C(M), we define the sets 

Ef = {lE M : f(x) = ll/IL}, £f = {ieM: f(x) = -||/|U} 

on which / attains its maximal modulus. Clearly, E^ U Ej ^ 0. Note that 
either E^ D E'jf = 0 or E^ D Ej = M, which happens if and only if / = 0. 

Proposition 9. Let M be a compact metric space, let f G C(M ) and let 
p G B(M). Then /iG J(f) if and only if 

IIMII = ll/lloo (3) 

and the Jordan decomposition of p satisfies 

p P {M \Ef) = 0 = p N (M \ Ef). (4) 

Proof. Let p G J(/). Then, clearly, (J3J) holds. Moreover, if 

Pp(M \ Ef) + Pn(M \ E^) > 0, 

then 

Kf) = [ f dpP ~ f fdpN 

Jp Jn 

< (p P (Ef) + p P (M \ Ef) + p N (E f) + p N (M \ E? )) ll/lloo 

= lll^lll ll/lloo = II/IIL 
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which contradicts //(/) = H/H^- Hence ([3]) holds. 

On the other hand, if ([2]) and (J3J hold, then 

M/) = [ f d hP ~ [ fd^N 

JEf JEf 

= (MEf) + M£f)) 11/11=0 = WII/lloo = ll/IIL 

so that /i G J(/). □ 

The following proposition provides an explicit formula for the semi-inner 
product on C(M). 

Proposition 10. Let M be a compact metric space and let f,g G C(M). 
Then 


(Ld)- 


II S'I |oo min{ min f(x) 
xe ep 


min —f(x)\ 

x€Ef 


with the convention min0 = oo. 

Note that = 0 = is impossible, and that therefore the right-hand 
side is finite. 


Proof. Since g is continuous, the sets E^ and E^ are nonempty and compact. 
Since f is continuous, it attains its minimum over Ef at some xf G Ef and 
its maximum over Ef 1 at some x^ G Ef 1 . As the Dirac measures d T p and 
S x n satisfy 5 x p G B(M) + and S x n G B(M) + , and because of 

6 x p(M\E^) = 0 = S x ,(M\E^) 

and ll&rf II = 114*11 = 1, Proposition 0 implies I^HoAf e J(g) and -|Mloo4* 
J(g). Therefore, Proposition [9] yields 


(f,g)- = inf(M/) : T e J(g)} < H^Hoo mm{8 x p(f), -4*(/)} 

= IMloo mm{f{x%),-f(x%)} = Halloo min{ min/(x),- max f(x)}. 

x£EP x€ Ef 

It is easy to see that no g G J(g) yields a lower value. □ 

When X = C( S^ 1 ) and A, B G /C c (R d ), explicit expressions for the sets 
Eff A _ aB and E^ A _ aB can be obtained using the following proposition about 
variational inequalities. 


Proposition 11. Let A E /C c (R d ), a* E A and x E R d . Then 

||x — a* || = dist(x, A) (x — a*, a — a*) < 0 for all a E A, (5) 

||a* — x|| = dist(A, x) (x — a*, a — a*) > 0 for all a E A. (6) 

Proof. Inequality flSJ) is standard (see e.g. [3 Proposition 7.4]), and (J6|) can 
be obtained by an analogous proof. □ 

We are now in the position to characterize the sets and Eff A _ ag . 

Proposition 12. Let A, B E /C c (R d ), and let <j a ,&b C £(R d ) be the corre¬ 
sponding support functions. 

a) If A = B, then E? A _ ffB = S"” 1 . 

b) If AC B, then Eff A _ ag = 0. 

c) Let A <f_ B. Then for any p E we have p E Eff A _ ag if and only if 

there exist a* E A and b* E B such that p = (a* — b*)/\\a* — b*\\ and 

||a* — 6*|| = dist(a*, B) = dist(v4, B) = dist# (A, B). (7) 

An analogous statement holds for the set Eff A _ ag . 

Proof. If A = B, then <j a = as, and hence 

E * A -a B = {P e : a A (p) ~ °b{p) = \Wa ~ *Lb 11oo} = 
which proves a). If A C B, then \\a A — ^bIIoo > 0 and a a ~ &b < 0, so that 

K a -* b = {P£ : a A (p) ~ a B (p) = \\<ta ~ HU} = 0, 

which is b). 

Let us show the equivalence c). Let p E Eff A _ ag . Using Proposition [2] we 
find 


dist h (A, B ) = 11 o 'a ~ ct b ||oo = <ta(p) ~ cr B (p) = sup (p, a) - sup(p, b) 

a£.A. b£B 

= sup inf (p, a — b) = sup inf cosZ(p, a — b)\\a — 6|| 

a&A b ^ B aeA b ^ B 

< sup cos Z(p, a — proj B (a))||a -proj B (a)||. 

aeA 
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By compactness of A and continuity of the above expression, there exists 
a* G A such that 

disttf (.A, B) < sup cos Z(p, a — proj B (a))||a — proj s (a)|| 

aeA 

= cos Z(p, a* — proj B (a*))||a* — proj B (a*)|| 

= cos Z(p, a* — proj B (a*)) dist(a*. B ) 

< cos Z(p, a* — proj B (a*)) dist(H, B) < dist(H, B ). 

Hence the above inequalities are, in fact, equalities, which enforces 

cos Z(p, a* - proj B (a*)) = 1, 

0 < distil(^4, B) = dist(H, B ) = dist(a*, B). 

Therefore, a* and b* : = proj B (a*) G B satisfy (J7|) and p = (a* — b*)/\\a* — b*\\. 

To show the opposite implication, let a* G A and b* G B satisfy (171) 
and set p = (a* — b*)/\\a* — b*\\. Note that dH) and the assumption A (jL B 
guarantee a* ^ b*. Using (jBj) and (J5]), we obtain 

(a* — b *, a) < (a* — b*, a*) for all a G A, 

(a* - b*, b) < (a* - b *, b*) for all 6 gB, 

so that 


sup (a* — b*, a) 

aeA 

sup (a* — b*,b) 

beB 


(a* -b*, a*), 
(a*-b*,b*}. 


Hence, using Proposition [21 we find 


<ta(p) - &b(p) = sup(p, a) - sup(p, b) 

a£A beB 

= ,| ,l b ,|, (sup (a* - b*, a) - sup(a* - b *, 6)) 

aeA b£B 

= ||a* -6*|| = dist h (A, B ) = \\a A - <t b \\oo, 


so that p G E* A _ aB 


□ 
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3 Existence and uniqueness of solutions 

In this section we apply standard existence and uniqueness results for the 
initial value problem 


x '(t) = x(0) = x 0 , (8) 

on a real Banach space X to the particular case of set differential equations 
in support function representation (121) . We first collect the necessary termi¬ 
nology and state a standard existence and uniqueness result for differential 
equations in Banach spaces from [6j. 

Definition 13. Let X be a Banach space, and let T>(X) be the family of 
all bounded subsets of X. The Kuratowski measure of non-compactness a : 
T>(X) —>■ R is defined by 

a(A ) = inf{d > 0 : A admits a finite covering by sets of diameter < d}. 

Definition [14] introduces standard classes of growth functions from [6]. 
The symbol D~ denotes the Dini derivative 

D~ pit) = liminf h~ 1 (p(t + h) — p(t )) 


of functions p : R — > R. 

Definition 14. We distinguish the following classes of growth functions. 

(UO) A continuous function u> : Rq —> Ro sa id to be of class Uq if the 
initial value problem 

p = u(p), p(0) = 0 

possesses only the trivial solution. 

(Ul) Let b > 0. A function uj : (0,6] x Mq —* M is said to be of class U\ if 
for each e > 0 there is a 6 > 0, a sequence L — * 0 + and a sequence of 
continuous functions pi : [t l} b] —> Mg" such that 

a ) Piifi) > dti for all i G N, 

b) 0 < pfit) < e for all i G N and t 6 (L, b], 
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c) there exists a sequence (5i)ig]N with 5* > 0 such that D pi(t ) > 
oj(t,pi(t )) + 5i for all i G N and t G (U,b\. 

The following existence and uniqueness theorem is an excerpt of [01 The¬ 
orem 4.1] applied in the present context. 

Theorem 15. Let ( X , || • ||x) be a Banach space, and let D C X, x 0 G D 
and r > 0 be such that D r = D n B r (x 0 ) is closed and convex. Let c > 0, let 
/ = [o. r] xf) r 4l be a continuous function satisfying 

\lf(t,x)\\x < C for all t G [0, T], x G D r , 

and let b := min{T, r/c}. Suppose that the subtangent condition 

fit, x ) G T D [x ) for all t G [0, b], x G dD fl B r (x 0 ) 

holds. Then the initial value problem (jHJ) has a solution <p : [0, b] —y D r 
provided one of the following additional conditions is satisfied: 

a) There exists a function u> : Rg —> Rg of class U 0 such that 

a(f([0,b] x A)) < u{a{A)) for all A C D r . 

b) There exists a function to : (0,6] x R + —x R + of class U\ such that 

{f{t,x) - f(t,y),x-y)- <w(t, \\x - y\\ x )\\x - y\\ x 
for all t G [0, 6] and x, y G D r . 

In case b), the solution is unique. 

When adapting Theorem [15] to set differential equations, we will fre¬ 
quently use the version 

f{t, a) G £(R d ) — Rg a for all t G [0, T], a G £(R d ), (9) 

of the subtangent condition to ensure that solutions do not leave the cone 
£(R d ) associated with /C c (R d ). 

Our first result is a Peano type theorem. 
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Theorem 16. Let f : [0,T] x £(R d ) — y E(R d ) be a continuous function, let 
A 0 £ /C c (R d ), and let r > 0. Then there exists c > 0 such that 

||/(^,cr)||oo < c for all t £ [0,T], a £ D r = E(R d ) fl B r (a Ao ). (10) 

Let b := min {T,r/c}. If, in addition, the subtangent condition ([9]) holds, 
then there exists a solution a : [0, b] —> D r of the set differential equation (J2]) 
in support function representation. 

Proof. Since balls defined in the maximum norm are always convex and 
£(R d ) is a convex cone, the intersection D r is convex. By Proposition |3l 
the set D r is compact, and the existence of some c > 0 such that (fTO holds 
is implied by the continuity of /. By Proposition [5] condition (|2]) implies 

f{t,cr) £ Pu(]R' i )( cr ) f° r all t e [0,6], a e D r . 

By compactness of D r and continuity of /, the image /([0,T] x D r ) is com¬ 
pact, and hence we have 

«(/([0, T] x A)) = 0 = a(A) for all A C D r , 

so that the compactness assumptions of Theorem fT5h ) are trivially satisfied 
with u)(p) = p of class U 0 . □ 

The next result is a Picard-Liudelof type statement. 

Theorem 17. Let A 0 £ /C c (R c/ ), and let f : [0,T] x S(R d ) —>■ E(R d ) be 
continuous and Lipschitz continuous in its second argument, i.e. we assume 
that there exists L > 0 such that 

\\f(t,<j A ) ~ f(t,cr B )||oo < L\\a A - ctbUoo = Ldist H {A,B) 

for all A, B £ /C c (R d ). If, in addition, f satisfies condition (0, then there 
exists a unique solution a : [0,T] —> E(R d ) of (J2]). 

Proof. As / is continuous and [0, T] is compact, we have 

K := SUp ||/(t,0A o )||oo < oo, 
te[o,T] 

and Lipschitz continuity of / yields 

c r := sup \\f(t,cr)\\oo<Lr + K. 

te[o,T], creBrCo-AolnsfiRd) 
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Because of 


- f(t,a),cr - a)_ = inf u{f(t,a)-f{t,a)) 

< JIMII \\f(t,a)-f(t, a) ||oo < L\\a - 

for all t G [0,T] and a, a G S(R a! ), and by the arguments in the preceding 
proof, all assumptions of Theorem fl5b) are verihed with r = 1, c = c\ and 
a;(i, s) = Ls, so there exists a unique solution cr 0 (-) : [0, b 0 ] —> S(R d )ni?i(cr J 4 0 ) 
of © with b 0 : = min{T, 1/(L + k )}. If 1/(L + k ) < T, the same argument 
yields a unique solution cr^-) : [b 0 ,b 0 + 6i] —>■ £(R d ) D Bi(a 0 (b 0 )) of the set 
differential equation with b\ := min{T — bo, 1/(2 L + k)}. 

Assume that b 0 + &i < T and that this construction can be repeated 
indefinitely with J2k=o bk < T for all JVeN. But then 

OO OO 

T >J2 bk = J2krr. = °°» 

fc=0 k =0 

which is a contradiction. Hence there exists a smallest index N G N such 
that &at = T. Concatenating the unique solutions ao,... ,<Tn yields a unique 
solution a : [0, T] —» £(R d ) of (jSD on the entire interval [0, T]. □ 

In contrast to the Picard-Lindelof type result above, the following state¬ 
ment fully exploits Theorem H5b) and the considerations from Section 12.41 
Roughly speaking, it states that uniqueness of the solution can be guaran¬ 
teed by controlling the relative velocity f(t, a a) —f{t, o~b) for two sets A, R G 
/C c (R d ) in only one critical direction that is given by a pair (a, 6) G A x B 
which realizes the Hausdorff distance of A and B. 

Theorem 18. Let f : [0,T] x E(R d ) —>■ ^(R 0 ') be continuous, let A 0 G 
/C c (R d ), and let r > 0. Let b, c > 0 and D r be as in Theorem\T^ let 

D' r := {A G /C c (R d ) : dist*(A, A 0 ) < r}, 

let u: : (0, T] x M + —y M be of class U\, and assume that the subtangent 
condition dHJ) holds. If, in addiditon, for any t G [0, b] and A, B G D' r with 
4 / B, there exist a & A and b G B such that p := (a — b)/\\a — 6|| is 
well-defined and one of the conditions 

||a — 6|| = dist(a, B) — dist(A, B) — dist#(A, B), 
fit , <ta)(p) - fit , 0b)(p) < w(t, dist#(A, B)) 
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Figure 1: Illustration of the geometric condition in Theorem [IH1 in an impor¬ 
tant special case. Let t G (0, T) and assume that there exist A, B G /C c (R d ) 
such that cr^ = f (t, a a) and = f(t,UB )■ Then the illustration depicts 
the relative velocity v r (p) = fit, cta)(p) ~ fit, cr){p) in the critical direction 
P=(a~ b)/\\a - b\\. 


and 

|| a — 5|| = dist (&, A) = dist(T>, A) = dist# (A, B), 
fit, <7 b)(-p) - f(t, <ta)(-p) < u(t, dist#(h4, B)) 
is satisfied, then there exists a unique solution a : [0,6] —> D r of fl2J). 

The geometric principle behind conditions (TTT1) and (fT2l) is depicted in 
FigureQ] for the case when f(t, a a), fit , o B ) G S. 


Proof. By Theorem [TIE we know that the desired solution exists. According 
to Theorem IT5b) . to ensure uniqueness, we need to verify that 

if it, (7 A) - fit, Cr B ), (7 a - a B )- < u(t, || (7a ~ Ob||oo) \Wa ~ 0b||oo 

for any t G (0, b] and a a, a B G D r . By Proposition fTOl this is true if and only 
if for any t G (0, b] and a a, (7b G D r , at least one of the inequalities 

min if it, a A ) ip) - fit , a B )(p)) < w(t, dist^(A, B )), 

P GE tT A -<T B 

min if it, a B )ip ) - fft, a A )(p)) < u(t, dist H (A, B )) 

PGE ‘ T B- rT A 
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is satisfied. If a a f <?b, this is, according to Proposition [[21 ensured by 
conditions m and dH} , which can be checked by addressing all possible 
relations A C B, B C A and A (jL B A B <f_ A between the sets A and B 
separately. If a a = erg, both inequalities are obviously valid. □ 

4 Hukuhara-type differentials 

In this section, we clarify that curves A : [0, T] —> /C c (R d ), which are sec¬ 
ond type Hukuhara differentiable, are time-reversed Hukuhara differentiable 
curves with the same derivative up to sign change. This insight has some 
important consequences. 

i) As Hukuhara differentiable curves can only grow in diameter, see Em 
Proposition 1.6.1], second type Hukuhara differentiable curves can only 
shrink in diameter, as claimed in the introduction. 

ii) As the support function representation of Hukuhara differentiable curves 
is Frechet differentiable, see [3j Lemma 4.1], this also holds for second 
type Hukuhara differentiable curves. Furthermore, by the same lemma, 
the Hukuhara and the second type Hukuhara differentials of a curve 
coincide with its Frechet differential (up to a sign change), whenever 
the Hukuhara type differentials exist. Therefore, set differential equa¬ 
tions based on both types of Hukuhara derivatives are special cases of 
the support function approach we presented. 

The notions of Hukuhara difference and Hukuhara differential are stan¬ 
dard. The concept of generalized or second type Hukuhara differentials goes 
back to [4j. Their use for set differential equations was investigated in [[131 ITT] , 

Definition 19. (Hukuhara differences and differentials) 

a) Let A,Be /C c (R d ). If there exists C E /C c (R d ) such that A = B + C, 
then C is called the Hukuhara difference between A and B, and we 
denote C = A Qh B. 

b) A curve A : [0, T] —$■ /C c (R d ) is called Hukuhara differentiable at t e 
(0, T ) with Hukuhara differential DnA(t ) G /C c (R d ) if the limits 

lim h~ l (A(t + h) Q H A(t)) , lim h~ l (A(i) Q H A(t — h )) 

h\0 h\0 

w.r.t. Hausdorff distance exist and equal DnA(t). 
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c) A curve A : [0, T] —* /C^R 0 ') is called second type Hukuhara differen¬ 
tiable at t e (0, T) with differential D* H A{t) e /C c (R a! ) if the limits 

lim(—/i) -1 (A(i) Q h Aft + h )), lim(— hy 1 (Aft - h ) Q H Aft)) 

w.r.t. Hausdorff distance exist and equal D* H A(t). 

The following proposition shows that second type Hnknhara differentiable 
curves are precisely those curves that are Hukuhara differentiable in the 
ordinary sense after time reversal. 

Proposition 20. Let A : [0, T] —> /C c (R d ) be a curve, and let B : [— T, 0] —> 
/C c (R d ) be given by Bft) = A(—t). Then A is second type Hukuhara differ¬ 
entiable att E (0, T ) if and only if B is Hukuhara differentiable at —t in the 
usual sense. In that case, the respective differentials satisfy 

D* H Aft) = —D H B(—t). 

Proof. The statement follows immediately from the identities 

lim(— h)~ l (Aft) Qh Aft + h)) = — lim h^ 1 (R(— t) Qh ( B(—t — h))), 

h\0 /i\0 

lim(—h) _1 (Aft — h) Q H Aft)) = — lim h^ 1 (i ?(—t + h) Qh B(—t)) 

h\0 /i\0 

for the Hausdorff limits. □ 


5 Example 

We conclude our paper with a simple, but instructive example, which il¬ 
lustrates that the usefulness of both types of Hukuhara derivative depends 
not only on the equation, but even on the initial value. Consider the set 
differential equation 


fr t °A{t) = cr Q - cr A (t), va(o)=va 0 (13) 

in /C c (R 2 ) with Q = [—1, l] 2 and A 0 = [ai,6i] x [a 2 ,& 2 ] C R 2 . The curve 


Aft) — e t A() + (1 — e t )Q 


(14) 
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Hukuhara differentiable solution corresponding Hukuhara differentials 




second type Hukuhara differentiable solution 



corresponding second type Hukuhara differentials 
1 

0 

-1 

-2 

-3 -2 -1 0 



not Hukuhara differentiable solution curve 



Figure 2: Solutions to set differential equation (1T3T) with three different initial 
values. The rectangles in the frames on the left are the values A(t), t = 
0, of the solutions. The rectangles in the top right frame are the 

Hukuhara differentials DnA(t), and the rectangles in the second frame on 
the right are the second type Hukuhara differentials D* H A(t ) at the same 
time points. The bottom right frame is empty, because the third solution 
curve is neither Hukuhara nor second type Hukuhara differentible. 


18 









































































Hukuhara differentiable solution 



corresponding differentials 




-2 
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0 2 4 6 



not Hukuhara differentiable solution curve corresponding differentials 




Figure 3: Solutions from Figure [2] and the corresponding differentials in 
support function representation. The differentials of the third curve cannot 
be interpreted as sets, but are well-defined as elements of S(R d ). 


19 














is a solution of (IT5j) . because 

lt a A(t) = e~ t (a q - o~a 0 ) = ctq - &A(t)- 

By Theorem |T7l the solution is unique. Clearly, the set Q is a globally 
asymptotically stable fixed point. 

By [3] Lemma 4.1], any Hukuhara differentiable solution of the reformu¬ 
lation 

DhA{P) — Q Qh A(t), A(0) = A 0 (15) 

of (fT3l) in set notation must coincide with this curve. Note that for many 
A G /C c (R 2 ), the right-hand side Q Qh A of (1T5li is not well-defined. Since 
<Tq — (Jao G £(R d ) if and only if 

max{&! — ai, b 2 — a 2 } < 2, (16) 

there does not exist a Hukuhara differentiable solution if this condition is 
violated. A computation shows that (1T6|) is sufficient for (ITT]) being a solution 
of first Hukuhara type. 

Proposition [201 however, shows that the curve (fTfl) can only be a second 
type Hukuhara solution, if ctq — <ja 0 G —S(R d ), which is equivalent with 

min {hi - cp, b 2 - a 2 } > 2, (17) 

and condition (TT7T) is sufficient for (fTTj) being a solution of second Hukuhara 
type. 

Figures [2] and [3] display solutions of (fT3]l with three different initial values 
Al = [2,3] x [1,2], Al = [0,3.5] x [-1.5, 2.5] and A 3 0 = [-1.5, 3.5] x [-0.5,0], 
Figure [2] depicts the sets as such on the left, ft is clearly visible that 
dist#(A(f), Q) —> 0 as t —> oo. The first curve is Hukuhara, but not sec¬ 
ond type Hukuhara differentiable, and the Hukuhara differentials are plotted 
in the top right subplot. The second curve is second type Hukuhara, but 
not Hukuhara differentiable, and the second type Hukuhara differentials are 
plotted in the middle of the right column. In both cases, the differentials 
converge to {0} when the state approaches Q. The third curve is neither 
Hukuhara nor second type Hukuhara differentiable, because it shrinks in the 
direction of the first and grows in the direction of the second axis. 

Figure [3] depicts the same three curves in support function representa¬ 
tion. The left column shows the evolution of the support functions, while 
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the right column shows the Frechet differentials along that curve. In this rep¬ 
resentation, the third curve can be treated as any other. The fact, that its 
differentials are elements of S(R d ) \ £(1R/ Z ) causes no problems. In all three 
cases, the derivatives converge to the zero function as the state approaches 
the equilibrium. 

We conclude that both types of Hukuhara differentiability only yield solu¬ 
tions for very special initial conditions, while the support function approach 
yields a solution that exhibits the expected behavior for any initial condition 
without any technical complications. 
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